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Abstract 



^ ' A semiclassical theory of linear response in finite Fermi systems, based 

on the Vlasov equation, and its applications to the study of isoscalar vibra- 
\Q ■ tions in heavy nuclei are reviewed. It is argued that the Vlasov equation 

can be used to study the response of small quantum systems like (heavy) 
in ■ nuclei in regimes for which the finite size of the system is more impor- 

tant than the collisions between constituents. This requires solving the 
linearized Vlasov equation for finite systems, however, in this case the 
problem of choosing appropriate boundary conditions for the fluctuations 
of the phase-space-density is non-trivial. Calculations of the isoscalar re- 
sponse functions performed by using different boundary conditions, corre- 
sponding to fixed and moving nuclear surface, are compared for different 
' multipoles and it is found that, in a sharp-surface model, the moving- 

surface boundary conditions give better agreement with experiment. The 
semiclassical strength functions given by this theory are strikingly similar 
to the results of analogous quantum calculations, in spite of the fact that 
, shell effects are not included in the theory, this happens because of a well 

■ known close relation between classical trajectories and shell structure. 

a i 

1 Introduction 



The Landau kinetic equation for Fermi liquids ^ contains some important 
differences when compared to the classical Boltzmann equation for a dilute gas, 
one of them is the presence of an effecive mean-field term. Thus, in Landau's 
approach, at least part of the force which is exerted on a particle by the other 
constituents of a many-body system can be approximated by an effective mean 
field. Another important feature of the Landau kinetic equation is the introduc- 
tion of an effective mass for quasiparticles, however here we neglect the difference 
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between bare and effective mass. Landau's kinetic equation allows also for a col- 
lision term which, nonetheless, in some cases may be neglected (see for example 
Sect. 4 of [3]). When this is done, we are dealing with a collisionless kinetic 
equation in the mean field approximation. Such an approximation had been 
considered long ago by Jeans (who gave the credit to Boltzmann) in the context 
of stellar dynamics @] and later by Vlasov in that of plasma physics [5] . Here 
we follow the use that has become common both in plasma and nuclear physics 
and refer to the collisionless kinetic equation in the mean field approximation 
as the Vlasov equation (see |H] for a discussion of historical priorities). 

Kirzhnitz and collaborators |7j extended the approach to non homogeneous 
systems and used it to study the possibility of collective excitations in the elec- 
tron cloud of heavy atoms. They pointed out that a main difficulty arising 
in finite systems concerns the boundary conditions to be imposed on the fluc- 
tuations of the phase-space-density. Moreover these authors also derived an 
interesting expression for the polarization propagator determining the linear 
response of these systems. Unfortunately the practical usefulness of this expres- 
sion is limited to rather special systems in which the constituents move along 
closed orbits. 

Another attempt to study the dynamic response of inhomogencous Fermi 
systems was limited to one-dimensional problems '8,- 

Bcrtsch P| argued that the Vlasov equation could be used as a starting point 
for a semiclassical theory of giant resonances in heavy nuclei. He pointed out 
that, in spite of being a classical equation of motion, this equation would not 
violate the Pauli principle, at least in a semiclassical sense. This is a consequence 
of the Liouvillc theorem. When applying this method to nuclei, however, one is 
faced with the problem of finite-size effects since a nucleon close to the Fermi 
surface is more likely to reach the nuclear surface than to suffer a violent collision 
with another nucleon. Therefore finite-size effects become more important than 
the collision integral and also in the case of nuclei it is reasonable to study the 
kinetic equation in the mean-field (or Vlasov) approximation, at least as a first 
step. 

Some remarkable progress on this problem has been made in the field of 
galactic dynamics: Polyacenko and Shukhman solved the linearized Vlasov 
equation for finite spherical systems in their study of the stability of collisionless 
stellar systems. In this context one of the main problems is that of determining 
a stable equilibrium distribution of particles (stars). Small deviations from the 
equilibrium distribution are characterized by eigenfrequencies that are purely 
imaginary in the case of unstable systems. Their approach has found several 
applications in the field ^1 EH • 

A similar solution of the linearized Vlasov equation for nuclear response has 
been derived independently in |14|. This solution turns out to agree with that 
of jlU) and gives a reasonable description of giant resonances in heavy nuclei 

E3 

Abrosimov, Di Toro and Strutinsky |16j used the same approach within 
a sharp-surface model in which the nuclear mean field is approximated by a 
square-well potential and used also different (moving-surface) boundary condi- 
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tions in order to extend the approach of |14| to low-energy surface vibrations in 
heavy nuclei. 

This paper is a review of work done in the last ten years, based on the 
approaches of Jl] and [T2|. In Sect. 2 both approaches are recalled, while in 
Sect. 3 several applications of the theory to the study of isoscalar vibrations of 
different multipolarity (monopole, dipole, quadrupole, octupole) in heavy nuclei 
are discussed. Finally, in Sect. 4, conclusions are drawn. The two Appendices 
contain some more technical material on the moving-surface response functions 
and on the Fourier coefficients that replace the quantum matrix elements in our 
semiclassical approach. 



2 Reminder of formalism 

2.1 Smooth surface 

In our semiclassical approach we assume that a (heavy) spherical nucleus in 
its ground state can be described by the following equilibrium phase-space- 
distribution 1 

/o(r, p) - J^j^HtF - h (r, p)) = F(h ) , (1) 
where d is the step function, tp is the Fermi energy, while 

Mr,p) = ^ + Vb(r) = e (2) 

is the quasiparticle energy and the equilibrium mean field Vo(r) is assumed to 
be spherical. In principle the mean field should be determined sclf-consistently 
as (Hartree approximation) 



V (r) = / drXr,r> (r'), (3) 
where v(r, r') is the effective interaction between quasiparticles and 

Oo(r) = / dp/ (r,p) (4) 



the equilibrium density of the nucleus (for sake of simplicity we do not take 
into account explicitly the spin and isospin degrees of freedom since this would 
only complicate the formalism without posing any new conceptual difficulty; 
the statistical weight 4 in Eq.Q accounts for these degrees of freedom). In 
practice we shall use instead a phcnomenological equilibrium mean field, which 
can be either of the Saxon- Woods shape or even a simple square- well potential 
of radius R = 1.2^3 fm. For the square-well approximation, several expressions 
can be evaluated analytically and this is one of the merits of our simplified ap- 
proach. Note that, contrary to the standard Fermi liquid theory, our equilibrium 
distribution depends also on the space coordinate r. 

x We use units h = c = 1 
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Next we assume that at time t = the system is subject to an external 
driving field of the kind 

V ext (r,t) = pS(t)Q(r). (5) 

Here (3 is a parameter determining the intensity of the external force, which 
is applied only for a very short time around t = 0, as described by the Dirac 
5-function, and Q(r) gives the space dependence of the external field. Typically 
we shall be interested in the multipole response of order L, for which 2 

Q(r) = r L Y LM (r). (6) 

The response of the system to an external force is described by the fluctuation 
of the phase-space density defined by 

/(r,p,t) = / (r,p)+5/(r,p,t) (7) 

or, equivalently, by its time Fourier transform 

/CO 
dte iut 5f(r,p,t) . (8) 
-oo 

Since 5f(r, p, t) vanishes for t < 0, we can suppose that uj has a vanishingly small 
imaginary part is to ensure the convergence of this integral when t — > +oo. 

The perturbed system is described by a time-dependent phase-space-density 
satisfying the mean-field (or Vlasov) kinetic equation 

^ + {/,M = 0, (9) 

where the braces are Poisson brackets. The time-dependent Hamiltonian h is 
given by 

h{v,p,t) = h (r,p)+5h{r,t), (10) 

with 

Sh{r, t) = V ext (r, t) + 5V mt (r, t) . (11) 

This expression shows explicitly that the extra force acting on a particle 
in the perturbed system has two components: one due to the external driving 
field V ext and an additional one due to the change in the interaction with the 
surrounding particles. This last term is given by 

5V int (r,t) = J dv'v{v,v') J dp'6f(r',p',t). (12) 

If the external force is sufficiently weak, the density fluctuation induced by 
it is small and we can consider only terms that are linear in [3. In this case the 
fluctuation 5f(r, p, t) satisfies the linearized Vlasov equation 

^f + { Sf,h } + {f ,6h} = 0, (13) 



2 For compression modes we'll be interested also in Q(r) = r L+2 YLM(r)- 
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Figure 1: Angular elements of the orbit. For a particle moving in a central force 
field, the angles a and are constant. 



or 

^f + {Sf,h } = -F'(h ){h ,Sh}. (14) 

From the mathematical point of view, if the (r, p) variables are used, this 
equation is a seven-dimensional differential equation (actually an integro-differential 
equation because of Eq. <|12|l ') containing partial derivatives with respect to time 
and to the six variables and pi. The time derivative can be eliminated by 
using the Fourier transform JSJ, while the properties of Poisson brackets suggest 
that some simplification might be achieved simply by making a change of vari- 
ables to generalized coordinates and momenta. The new coordinates should be 
chosen in such a way to include the maximum number of constants of motion of 
the unperturbed Hamiltonian Hq so that the corresponding variable will not con- 
tribute to the Poisson bracket. For motion in a central force field a convenient 
set of generalized coordinates is (e, A, r, a,0, 7), where e is the particle energy 
(|2Jl, A the magnitude of its angular momentum, r the radial coordinate, and 
a,f3,j the Euler angles associated with the rotation of the frame of Cartesian 
coordinates necessary to align the z-axis of the lab frame to the particle angular 
momentum A and the y-axis of the lab frame with the r vector specifying the 
instantaneous position of a particle with respect to the force centre (see Fig. 1). 

Since four of the six new coordinates are contants of the motion, Eq. I|14|l 
simplifies considerably and becomes 

dSf d Sf . d8f . . t „d5h . dSh., 
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By using 



7 = —, (16) 

TJIT 

f = ±v r (e,X,r), (17) 

with 

^,A,r) = ^[e- W -^] (18) 
the magnitude of the radial velocity, we obtain 

-iu>6f±v r (e,\,r)^- + J = F (e) [±v r (e, A, r) — \- 19 

or mr z 07 or 07 

for the Fourier-transformed linearized Vlasov equation. 

This equation still contains partial derivatives with respect to the two time- 
dependent variables r and 7, but the derivative with respect to 7 can be elim- 
inated by means of an appropriate partial-wave expansion. The usual partial- 
wave expansion 

5f(r,p,u) = '%25f L M(r,p,u)Y L M(4,<p), (20) 

LM 

where (1?, tp) are the polar angles of the vector r, is a first step in this direction, 
but it does not solve the problem since the coefficients <5/LM( r , P, w ) still depend 
on the direction of the vector p, however the following well known transforma- 
tion property of the spherical harmonics Ylm{"&, f) under the rotation specified 
by the Euler angles a,/3, 7 can help us (see e.g. (1.7,, p. 28) 3 : 

L 

Y LM (#,v) = (^W^Vijv^V). (21) 

N=-L 

In the new reference frame the particle is on the y axis, so i?' = ? and 
tp' = \ and the only time-dependent angle on the righ-hand-side of Eq. (|21|l is 
7. The functions T> MN (af3,j) are the coefficients of the rotation matrices and 
their explicit 7 dependence, of the kind e~ lN ^ , can be exploited to eliminate 
the 7-derivative in Eq. I|19|l , giving the following one-dimensional equations 



0_ 

Or 

with 



Z^MAT T A N S fMN ~ Bajn 1 (22) 



a 1 \ * w A , . 

A N (e, X,r,Lu) = - (23) 

u r (e, A, r) v r [e, A, r) mr z 



and 



B^ N (e, \,r, U ) = F'(e)[^ ± -£^-^][PQ LM {r) + SVffi{r, U )] , (24) 



3 In this rather technical aspect the derivation of 1141 differs by that of II 1)1 . which is less 
straightforward (cf. also llll l. 
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for the coefficients Sf^^(e, A, r, u>) of the expansion 4 

5/(r,p,«)= ^ [5/ A L f + (6,A,r > a;)ff(p P ) + */^(e,A J r,w)e(-jj r )] 

LAIN 

(v L MN {a(] n )yY LN {~). (25) 

In Eq. H24|) the functions QLM{ r ) an d SVjJ^^r^uf) are the coefficients of a 
multipole expansion similar to (|2U[I for the external driving field and for the 
induced mean-field fluctuation, respectively. In Eq. (|25() instead, the ^-functions 
are the usual step function 9{x) = lfor x > and 0(x) = otherwise, while p r 
is the radial component of the particle momentum: p r ~ ±mv r - 

Thus, by making an appropriate change of variables (and by taking the 
time Fourier transform), the initial seven-dimensional differential equation l|14fl 
has been reduced to the system of two (coupled) one-dimensional differential 
equations l|22|) . These two equations, involving the distributions Sf + and 8f~ 
of particles with both signs of the radial velocity, are coupled by the mean-field 
fluctuation, as shown explicitly by the following expression: 

!«/£» ( e . A . >•» + «/fev(e. \r',u)\. (26) 

For most effective interactions vl(t, r'), because of this coupling the solution of 
Ea. ()22j) can only be given in implicit form, however, an explicit solution can be 
obtained if we neglect the term SV™ K \{r^uj) in Eq.J2U). Following Ref. ^3], we 
refer to this as the zero-order approximation and recall here the solution details. 

In order to solve Eq. (|22f> . we must first specify the boundary conditions 
satisfied by Sf + (r) and Sf~(r) at the turning points n and ?- 2 . The boundary 
conditions used in |14j were: 

^/m7v( £ ' \ r i,v) = £/MAr( e > A > r i> CJ )' ( 27 ) 
(5 /Miv( e ' \r 2 ,u) = Sf^f N {e,X,r 2 ,uj) , (28) 

their physical meaning is that, at the turning points the radial motion of the 
particles simply reverses. For the square-well potential, the condition at the 
outer turning point implies a mirror reflection of the nucleons on the static 
equilibrium nuclear surface. As we shall see, in the sharp-surface case there 
are reasons for modifying this boundary condition, however for the moment we 
assume a diffused surface and determine the solution of Ea. (|22|l by using the 
boundary conditions (|27|l and (|28(l . The solution can be written as (in slightly 
simplified notation) 

*/j&(e,A,r,w) =e ±i ^^[[\r'B^ ± N (r')e^^ r '^ +C±(e,A,w)], (29) 

±ri 

4 Note that, contrary to what is done in Ref. 1141 . here the factor ij,iv(^: 5) i s n °t included 
in our definition of the coefficients Sf^ N . 



7 



with r 

4> N {r,u) = -ij dyA N (y) = wT(r) - N-y(r) , (30) 

r{r) -iy^yy (31) 

The functions C±(e, A, lu) play the role of integration contants and are deter- 
mined by the boundary conditions. The inner boundary condition l|27|) implies 

C- = C + , (33) 

while the outer boundary condition (|28() implies 

/ 2 dr'B^ N {r')e- t ' t ' N{ - r '^ ) + C+] = 

Jri 

-i<Mr 2 , w )[ /" 2 dr'S^^Oe^^''^ + C+] . (34) 



e 



Defining 

gives 
that is 



D±(e,\,u>)= f ~ dr'B L ± N {r')e 

Jri 



(35) 



e 2 



2i0jv(r-2,w) n n 

C+(e,A,w) = ± , ~ =C_(e i A,cj). (37) 

The most interesting property of the solution l|29|) is its pole structure in 
the complex-w plane, which is entirely determined by the poles of the functions 
C±(e, A, u>), that is, by the vanishing of the denominator in Eq. <|37[) . This 
happens whenever 20jv(^2,w) = n2n, with integer n. This is the point where 
the finite size of the system plays a crucial role since the eigenfrequencies of 
the density fluctuations in the zero-order approximation are determined by the 
condition 

w[2r(ra)] - N[2j(r 2 )] = n2ir. (38) 

and the period of radial motion T(e, A) = 2t(j2) depends on the size of the sys- 
tem. For motion in a central potential, these eigenfrequencies are characterized 
by the two integers n and N: 

Lo nN (e, A) = nujQ + Nuo v , (39) 
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with 



2-7T 

"o(e,A) = ^, (40) 

Me,A) = ^M. (41) 

the frequencies of radial and angular motion for a particle with energy e and 
angular momentum magnitude A. The eigenfrequencies <|39|) can be related to 
the difference of single-particle energy levels e n i in a central potential since, for 
large quantum numbers (UE|, p. 579), 

en'V - e»«r « (n' - n")-^ + (/' - f')-^ = 7iav-n'> . (42) 
an at 

Thus the integer n can be interpreted as the difference between radial quan- 
tum numbers in a single-particle excitation and the integer N as the difference 
between the corresponding orbital quantum numbers. 

The time-dependent density fluctuations can be obtained from the functions 
8fi^r(e,X,r,Lo) by contour integration in the complex-w plane. For the zero- 
order fluctuations we obtain: 

Sf^(e,X,r,t) =0 for K0, (43) 

1 f + °° 

= ^J due^Sf^^X^u) for t > , (44) 

with 

rj?0L±/ \ \ arplf \ ±id> nN (r) QnN^,X) . . 

°Jmn ( e A?",w) = -0F (e) > w niv e > — , (45) 

' — ' cj — W„tv + 

n=— oo 

nJV (r) = w nJV T(r) - JV 7 (r) (46) 

and 

Q nJ v(e, A) = — !-r f ° dr Q ^ M ^\ cos[0 nJV (r)] . (47) 

This result is obtained from Eq. I|29[l by closing the integration path in the 
lower part of the complex-w plane in the integral (|44|) . 

The coefficients (|47|) correspond to the classical limit of the quantum matrix 
elements of the external field (JSJ) T9*|. 

Thus we have seen that, for a spherical nucleus, the linearized Vlasov equa- 
tion can be solved explicitly in the approximation in which the mean-field fluc- 
tuation is neglected. This zero-order solution can be used as a starting point for 
solving also the more general problem in which the mean field fluctuation 5V mt 
is taken into account. The zero-order solution is most conveniently expressed in 
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terms of a semiclassical propagator ( obtained from Eq. 1451) . see ^31) which is 
analogous to the quantum particle-hole propagator 

s 2 +°° L f f 

£°(r,r» = ^— £ E \Yln(~)\ 2 / deF'(e) / dXX 

n=-ooN=-L ^ ^ 

1 COs[(/> w jy(r)] UJ nN cos[0„at(/)] 

T(e, A) r 2 v r (e, A, r) w„jv - (ui + is) r' 2 v r (e, A, r') ' 

Taking into account also the mean-field fluctuation (|12fl gives the collective 
response of the system and leads in general to an integral equation which is 
analogous to the RPA equation for the quantum propagators 

D L (r,r',u) =D° L (r,r',uj)+ / dxx 2 / dyy 2 D° L (r,x,uj)v L (x 1 y)D L (y,r l ' ,u) . 



(49) 

The only differences between the present expression and the quantum result 
are that here the propagator D° is given in semiclassical approximation and, 
of course, the exchange (Fock) term is missing. The integral equation (|49(l can 
be easily solved numerically (the zero-order propagator is actually simpler than 
suggested by Eq. 148|) since, for small values of u>, the infinite sum over n can 
be approximated with a sum over a few terms around n — with very good 
accuracy). The collective (multipole) response function is then given by 

TZ L (u)= [ drr 2 [ dr 'r' 2 Q LM (r)D L (r,r' \w)Q LM (r') , (50) 



with Dl solution of Q49[) . and the corresponding strength function by 

S L {Tluj) = - -IxoR-l (Uuj) (51) 

(for a spherical system the response is independent of M) . 

However there is a special case in which also the collective solution of the 
linearized Vlasov equation can be obtained explicitly. This happens if the inter- 
action between particles is supposed to be of the separable multipole- multipole 
type: 

v(r, r') = E ^r L r' L Y LM (r)Yl M (r') (52) 

LM 

and the external field is also of the multipole type 10. In this case Eqs. (|49|l 
and (|50|l give immediately 



KlH= L 'J, r , (53) 



with the zero-order response function TZP L (io) given by 



dXX 



v r {e, A, r) 
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Moreover the fluctuation Sf can also be obtained in explicit form, like Sf°. For 
the separable interaction l]52[l. Eq. I|26|) for the mean-field fluctuation gives 

SVi%(r,u,) = /3K L r L TZ L (u>) , (55) 

while Eq. (|2~lf) gives 

(e, A,r, W ) = F'(e)[|; ± -^L-^}^ + 0K L r L K L (w)} , (56) 

and, from Eq.C3, we get 5f*f N (e, A, r, w) = 6f°£±(e, A, r, w)[l + k l K l (uj)] or, 
by using Eq. (J53J, 

^( e ,A,r,a,) = ^ ± ^^ ) - (^) 



2.2 Action-angle variables 

Up to now we have assumed a spherically symmetric equilibrium mean field, 
however the method outlined here for the solution of the linearized Vlasov equa- 
tion is valid also for a wider class of physical systems. This method, which is 
based on the use of generalized coordinate and momenta in order to simplify the 
Vlasov equation, can actually be used for all systems which are described by an 
integrable equilibrium Hamiltonian /i (r,p). Such a Hamiltonian includes also 
some deformed systems (201 • For integrable systems, it is convenient to intro- 
duce action-angle variables (I, «&) instead of (r, p) since in this case the action 
variables I a are constants of the motion, while the angle variables <& Q are linear 
functions of time (see for example Ref. j2J, p. 457). An important property of 
these variables is that the motion is periodic in the angle variables with period 
2tt. Consequently the field felt by a particle that is moving along a trajectory 
determined by an integrable Hamiltonian can be Fourier expanded as 

^ t (r, W )=/3^Q n (I)e in -* ) (58) 

n 

and 

SV int (r, «) = ^ 5V n {I, Lu)e ln -* (59) 

n 

where n is a three-dimensional vector with integer components. 

The phase-space-density fluctuation can also be expanded in the same way: 

<J/(r,p,u,)=^/ n (I, W )e in -* (60) 

n 

and the linearized Vlasov equation gives the following equation for the coef- 
ficients <5/ n (I, uo): 

5f a (I, u) = F'(e) \pQ u Q) + SV**(I, n)l _ " ^ . , , (61) 
L J n • lu — (lu + is) 
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where the vector cD has components 

_ dh (I) 



(62) 



Again, Eq. (|61|) gives only an implicit solution of the Vlasov equation since 
the mean-field fluctuation SV lnt depends on 6f. If the term 5V mt is neglected, 
then Eq. (|61() is an explicit solution that, in the case of spherical systems agrees 
with the zero-order solution (|4"5|l . Generally speaking, the vector uj has three 
components, but in spherical systems there are only two basic eigenfrequencies: 
wo and u v . This is because spherical systems are over-integrable and this implies 
that one of the angle variables is also a constant of the motion. The coefficients 
Q n (I) in Eq. (|6*TJl are given by 

Qn(I) = 5^3 / d*e-*-*Q(T) (63) 
and correspond to the quantum matrix elements of the operator Q(r). 



2.3 Sharp and moving surface 

With the aim of establishing a link between the present microscopic theory of 
nuclear excitations and the macroscopic description given by the liquid-drop 
model (see e. g. [IB], Appendix 6A) the authors of ^E] studied in detail the 
model in which the equilibrium mean field is approximated by a square-well 
potential: Vo(r) = — Vb#(i? — r). They noticed that in this case the boundary 
condition (|28|l corresponds to a mirror reflection of nucleons when they reach 
the static nuclear surface at r = R. With that boundary condition the sharp- 
surface model allows only for compressional excitations, while a liquid drop has 
both surface and compression modes. They argued that, in order to allow also 
for a microcopic description of surface modes, in the sharp-surface case, the 
boundary condition should be modified and they proposed to replace it 
with the following (moving-surface) boundary condition: 

Sf M + N (e, A, R, u) = 6fjf N (e, A, R, u) + F'(e)2mv r (e, A, R)iu5R lm (lo) . (64) 

Whenever the same symbol is used, we put a tilde over the quantities evaluated 
with the moving-surface boundary condition 164|) . to distingush them from the 
corresponding quantities satisfying the fixed-surface boundary condition 
The physical picture behind this new boundary condition still corresponds to 
a mirror reflection of nucleons at the nuclear surface r = R, but in a reference 
frame that is moving in the radial direction at a speed R^, ip, t). The radial 
momentum transfered by the moving surface to the impinging nucleon will differ 
from that occurring in a collision with the static surface and this modifies the 
momentum distribution of nucleons according to (|64|l . in first approximation. 
Formally, in this approach the nuclear surface is allowed to vibrate (for t > 0) 
according to the usual liquid-drop model expression: 

JZ(0, <p,t)=R+Yj $RLM(t)Y LM ($, ip) , (65) 

LM 
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giving _^ 

R(#,cp,t) = ^6R LM {t)Y LM {fi,<p) ( 66 ) 

LM 

for the surface speed and 

R(0,<p,w) =J2- iu > SR LM(u)Y L M($,v) (67) 

LA/ 

for its time Fourier transform, thus leading to the extra term in Eq. (1641) . 

The new collective coordinates 8Rlm(w) are still to be determined. In JS] 
this has been done by recalling that in the liquid-drop model a change in the 
curvature radius of the surface results in a change of pressure given by (see Eq. 
(6A-57) of HH|) 

SV(R, 9, V ,u) = Y, Cl SRl *£ u) Y lm (9, y) , (68) 

LM U 

with the restoring force parameter Cl that can be related to the phenomenolog- 
ically determined surface tension parameter a ~ 1 MeVfm~ 2 . If the Coulomb 
repulsion between protons is neglected, this relation is simply 

C L = <jR 2 (L~1)(L + 2), (69) 

while taking into account also the Coulomb interaction gives an additional con- 
tribution to Cl (see [T2|, p-660). The pressure fluctuation (|6~5|) can also be 
related to the appropriate component of the pressure tensor (generalized to 
Fermi liquids, see HU, Sect. 74) 

6V(R,e,tp,cj) = J dpmv 2 (8f(R lP ,uj) - F'(e)SV int (R,u)) . (70) 

By equating the pressure fluctuations given by Eq. (|68|l and (|70|l , the collective 
coordinates 8Rlm{^>) can be related to the phase-space-density fluctuation: 

8ir 2 R 2 1 



5R LM (w) = — \Yln{-x,-x) de d\Xmv r (e, A, R) 

L N=-L J J 

5f^\{t,\ 1 R,u) + 8flf N {e,\R,u)~2F\e)5Vlt{R 1 u)\ . (71) 



The internal part of the mean-field fluctuation 8Vj^(r,ui) in the moving-surface 
case will be specified better in Appendix A. This can be done most easily by as- 
suming a separable interaction of the kind 1)52(1 . We have already seen that such 
an effective interaction leads to simple analytical expressions for the multipole 
response function and for the solution of the linearized Vlasov equation with 
fixed-surface boundary conditions. The same happens also in the moving-surface 
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case, although the final expressions are somewhat more involved. Since the ex- 
plicit derivation of the moving-surface multipole response function is rather 
lengthy, we report here only the final result 

n L {uj) = n L {uo) (72) 
i [xiM + £A« L i^ft°H] 2 

1 - n L nl(io) [C L - xl{u>)][1 - k l TZ° l (a;)] + K L b&(u>) + -^AR L } 2 ' 

and outline it in Appendix A. The response function 1Zl(uS) in the equation 
above is still that given by Eq. (|53|l . while the functions xl{u) and Xm^) are 
defined as 5 

Sir 2 



N=-L 

de [ d\\2F'(e)cot[<f> N (R,cj)][mv r (e,\,R)] 2 u. (73) 



and 



i o 2 ^ 
n / , 1 B7T v — > I , 7T 7T , 



N=-L 



de 



dXXmv r (e, A, i?) [<5/°f +(e, A, i?, w) + «/j&7(e, A, i?, u)] . (74) 



3 Isoscalar excitations in heavy nuclei 

In the remaining part of this paper we review some applications of the theory 
outlined in the first part to the study of isoscalar nuclear response . Our approx- 
imation for the mean field ( a square-well potential) is not particularly realistic 
since it neglects the surface diffusion and the assumed residual interaction of 
the multipole-multipole type is also rather special, however these approxima- 
tions have the advantage of leading to simple analytical formulae for the nuclear 
response functions. Our semiclassical approach shows that some feature of the 
nuclear response that are usually ascribed to quantum effects can be understood 
in terms of classical concepts like nucleon trajectories. The correspondence be- 
tween shell effects (not included in our treatment) and the properties of classical 
trajectories helps to shed a new light on the nuclear response at low energy. 

Our starting point is the zero-order response function 1|54[) . By using the 
explicit expression (|45|l of the phase-space-density fluctuations, with F'(e) — 
4 r 5(e.F — e), this equation gives 



L +oo 



N— — L 7i— — oo u v J 

(75) 



3 In Refs. [23 24 , a different normalization of has been used. 
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with x = sina, cosa = X/(pfR), s — u>/(vp/R), vf = Pf/™, 

nn + Na(x) 

SnN{x) = , (76) 

Qnkx) = lf dr r \ cos[0 nJV (r)] , (77) 

and 

Instead of the frequency u>, we have introduced the dimensionless parameter s 
and instead of the particle angular momentum A, the parameter x = sina. 

The Fourier coefficients (|77|) can be easily evaluated explicitly, the expres- 
sions needed here are grouped together in Appendix B. 

In terms of the new dimensionless variables, the auxiliary functions Xl an( i 
Xl appearing in eq. I|72[) read 



nN J ° 

and 



Xl(s) = J^R L J2 C LN C dxX*S nN { X ) "- >^*";\ (7!). 

47r J s + ie-s nN (x) 



X L (s) = -^'E^J^ } . (80) 

nN 

In the last equation we have used the pole expansion of the cotangent: 



cot Z = 



Z 1 y — 7 



z - 

n— — oo 



3.1 Monopole response 

This channel corresponds to a compression mode. For these modes the radial 
dependence of the external field is not given by Eq. |JB}, but rather by 

Q(r) = r L+2 Y LM {v). (81) 

As a consequence the zero-order response function (|75|1 willl involve the 
Fourier coefficients QnN^( x ) instead of Ql^ N (x) and the same is true for the 
auxiliary function Xl( s ) m ^1- GHl- Accordingly, the moving-surface response 
function (|72|l is also slightly changed to 

K L (w)=R L (u) (82) 
1 [xl^) + ^AK L R L ^Ul{u)f 
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i - K L nl{u>) \c L - X lW - K L n° L (u)] + K L [xl(") + i; ARL+ 

For L — this expression gives the collective moving-surface strength func- 
tion shown in Fig. 2 (full curve, E = hui). 
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E (MeV) 



Figure 2: Monopole strength function in three different approximations: zero- 
order (dotted curve), collective with finite value of Kl (solid) and collective with 
kl = (dashed). Note the change in the vertical scale. 



The collective strength function of Fig. 2 has been calculated by assuming 
a residual interaction 

VL=o(r,r') = K L=0 r 2 r' 2 , (83) 

with a value of kl=o = — 2 x 10~ 4 MeVfm -4 . This parameter has been deter- 
mined by fitting the experimental position of the giant monopole resonance in 

208p b 

The dotted curve instead shows the zero-order strength function (propor- 
tional to the imaginary part of the response function 1Z° l (uj), which is similar to 
the quantum single-particle response function). 

Finally, the dashed curve shows the collective moving-surface response given 
by Eq. I|82|) with L = and kl — 0. If kl = 0, the frequency of the collective 
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monopole vibration is determined by the solution of the equation 

C l -Xl(p) = 0- (84) 

It has already been pointed out in that this approximation gives a very rea- 
sonable description of the position (including the ^-dependence) of the isoscalar 
giant monopole resonance in heavy nuclei. 

Another interesting feature of the monopole response pointed out in is 
that the zero-order strength function vanishes for u> < ttvf/R- As a consequence 
of this fact, within this model, there is no Landau damping of the collective 
monopole mode. This absence of Landau damping is in qualitative agreement 
with the results of analogous quantum calculations |2s3 OSl- The very small 
width appearing in Fig. 2 is due to our use of a fnite value of the infinitesimal 
parameter e ( for numerical reasons, we have used e = 0.1 MeV). 

We have checked numerically that the collective state shown in Fig. 2 ex- 
hausts about 99% of the monopole energy-weighted sum rule, which is given by 

dE ES(E) = -^-—AR 2 . (85) 
l(j7r m 



3.2 Dipole response (translation and compression modes) 

It is well known that the mean-field approximation violates the translation in- 
variance of the nuclear Hamiltonian and that this results in the appearing of 
spurious strength in the isoscalar dipole response. Hence the isoscalar dipole 
channel, excited by the external field with L — 1, is usually not interesting 
because it should correspond to a simple translation of an unexcited nucleus, 
while in the mean-field approximation this pure translation is replaced by a 
spurious excitation of the nuclear centre of mass bound by an unphysical force. 
However the corresponding compression mode, excited by the field 

Q(r) = r 3 F 1M (r), (86) 

has received considerable attention because of the possibility of obtaining from 
it additional information about the compressibility of nuclear matter [2Hj • Since 
the external field l|86(l can excite also the centre of mass, the problem of subtract- 
ing the unwanted spurious strength from the corresponding response function 
has usually been dealt with by using an ingenious trick due to Van Giai and 
Sagawa |2H]. These authors suggested that, instead of studying the response to 
the external field (j86|l . one should look at the response to an effective external 
field of the kind 

Qeff(r) = (r 3 -vr)Y 1M (r), (87) 

where r\ is a parameter determined by the condition that, under the action of 
the external force, the centre of mass should remain at rest. 

The moving-surface theory of allows for a different approach to the 
problem of evaluating the intrinsic response associated with the field (|8fi|> . It 
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is clear from the expression (|69l) of the restoring-force parameter Cl that this 
parameter vanishes for L = 1. This means that in this case there is no restoring 
force, hence the moving-surface boundary condition seems to be able to readjust 
the translation symmetry that is broken by the mean-field approximation. This 
statement can be easily verified by looking at the form taken by the response 
function i|82[l when the residual interaction is neglected. If we put kl = in 
Eq. JH2J, the moving-surface response function becomes 

^ =1 ( W )=^H-M^, (88) 

or (cf. Eq. (3.1) of Ref. 0) 

which has the behavior expected for a free particle. Since this response function 
has no poles for u ^ 0, it does not give spurious dissipation at positive ui. 

The translation invariance of the model when the residual interaction is taken 
into account is less obvious, however it has been shown in by using sum 
rule arguments, that no spurious strength is added to the zero-order intrinsic 
response even when a residual dipole-dipole interaction is included. Thus we are 
confident that the present semiclassical theory with moving-surface boundary 
conditions correctly separates the intrinsic from the centre of mass excitations 
and do not need to use the effective operator (|57Jl . 

For sake of simplicity, here we discuss in some detail only the zero-order 
response, in which the residual interaction is neglected. This approximation 
corresponds to treating the nucleus as a gas of non-interacting fermions confined 
to a spherical cavity with perfectly reflecting walls that are allowed to translate 
freely. The residual interaction changes the compressibility of this nuclear fluid 
and its effects on the response have been evaluated in Ref. [231 - 

We first need to give a slight generalization of the zero-order response func- 
tion l|75fl . by defining the functions 

9A R j+k A 2 +^ /-i 2 ,JQW(x)/Ri)(QW(x)/R'<) 



*U« = ^— E °In E f d ™ 

16tt e F J~-> r Jr^^h 



s u n{x)- 



N=-L 



s - SnN(x) + ie 

(90) 

In terms of these new functions, the fixed-surface zero-order response to the 
external field i|86|) is determined by the function TZ\ ^{s) : while TZ® 13 (s) has a 
direct physical interpretation as the displacement of the nuclear centre of mass 
induced by the external field 23 . Similarly 1Z® n (s) gives the response to 
the field rY 1M (r). 

It has been shown in |23| that the (zero-order) moving-surface response func- 
tion for the field H8(jfl can be written as 

£l,33(*)=£° ro .(*)+£Lr(«). (91) 
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Figure 3: Comparison of our energy- weighted strength function with data from 
Ref. j.'iOj . The curve shows the strength function in the case of vanishing residual 
interaction, i. e. for a confined Fermi gas with an incompressibility of K = 200 
MeV. 



with 



and 



7?° 

' vntr 



V I 1,33 V, J 4?r 2 s 2 \ 



1 - 



i S 2 
2 A 



7? 1 



ju 1 



1 „2 Miiilfl 
~ S All 



The moments 



/ ^[--Im^? jfc ( S )] 
Jo n 



(92) 



(93) 



(94) 



are defined in terms of the generalized fixed-surface response functions H9Ufl and 
they can be easily evaluated, giving M\ x = f^r^ and M\ 3 = i? 2 A^^. An 
essential property of the intrinsic response function (|93|l is that its limit for 
s — > is finite, so it has no pole in u> = 0. 

A properly normalized and energy-weighted strength function associated 
with the intrinsic response function 1)93(1 is shown in Fig. 3. It is interesting to 
note that the simple model used here qualitatively reproduces the experimental 
data of 30 , in particular the double-peak structure of the dipole compression 
mode. Clearly at this level we can only hope in a qualitative agreement since 
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the present model has several unrealistic aspects (sharp surface, no residual 
interaction, etc.). 

3.3 Quadrupole response 



A=208 
Quadrupole 





Figure 4: Quadrupole strength function for a hypothetical nucleus of A = 208 
nucleons. The solid curve shows the moving-surface response, while the dashed 
curve gives the fixed-surface response. The dotted curve shows the response in 
the zero-order approximation. 



Figure 4 shows the quadrupole response function given by Eq. (|72|l (solid 
curve) with the strength of the residual quadrupole-quadrupole interaction de- 
termined by jTTT) 

k l =2 = -1.0 x 10~ 3 MeVfm" 4 . (95) 

The value of this parameter has been fixed by the requirement that the 
position of the giant quadrupole resonance (GQR) in our hypothetical nucleus 
of A = 208 nucleons agrees with the experimental position of the GQR in 208 Pb. 
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The obtained value turns out to be about twice that given by the Bohr-Mottelson 
prescription (|T^|, p. 509). Taking into account the fact that our equilibrium 
mean field has a different shape (square well rather than harmonic oscillator) 
and that we are assuming a semiclassical framework, this kind of agreement 
looks reasonable. While the position of the GQR is well reproduced by the 
appropriate value of /?2, its width is severely underestimated by our theory. 
This is a well known limit of all mean-field calculations in which the width 
is generated only by Landau damping, including a collision term in our kinetic 
equation would increase the width of this resonance. It is interesting to compare 
the moving surface response with the fixed-surface one (dashed curve in Fig. 
4): contrary to the fixed-surface response function, the moving surface response 
displays a low-energy bump whose exact position is determined by the value of 
the surface tension parameter a of Eq. (|69l) . Thus the moving-surface theory 
does reproduce both systematic features of the quadrupole nuclear response 
that are the GQR and the low-energy surface excitations. Quadrupole response 
functions calculated for other values of A, corresponding to other medium-heavy 
spherical nuclei, are qualitatively similar to the case shown in Fig. 4. 
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3.4 Octupole response 

For L = 3, Eq. i(72|) gives the moving-surface collective octupole response func- 
tion. A detailed study of this case has been made in Ref. It is interesting 
to look first at the fixed-surface zero-order response. It can be seen from Fig. 5 




E (MeV) 



Figure 5: Semiclassical octupole strength function analogous to quantum single- 
particle stregth function. Calculations are for A = 208 nucleons in a square-well 
potentialof radius R = 1.2A 1 / 3 fm. 

that, for A = 208, the single-particle octupole strength is concentrated in two 
regions around 8 and 24 MeV. As pointed out already in |14|. in this respect 
our semiclassical response is strikingly similar to the quantum response, which 
is concentrated in the lhu and 3hui regions. This concentration of strength is 
quite remarkable because our static distribution, which is given by Eq. Q , does 
not include any shell effect, however, because of the close connection between 
shell structure and classical trajectories expressed by Eq. (|42|l . we still obtain 
a strength distribution that is very similar to the quantum one. 

The collective moving-surface octupole strength function is shown instead 
in Fig. 6 (solid curve). 

Again we obtain a qualitative agreement with experiment and with the result 
of analogous quantum calculations. Like for the quadrupole case, agreement 
with experiment can be obtained with a residual interaction parameter about 
twice that given by the Bohr-Mottelson prescription. The rather broad double 
hump on the low-enegy side has been interpreted as a superposition of surface 
vibrations and of the low-enegy component of the giant octupole resonance. 
Within the present semiclassical theory, it can be shown that, for L = 3, the 
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Figure 6: The solid curve shows the octupole strength function given by the 
moving-surface solution (|72l) . while the dashed curve gives the corresponding 
fixed-surface response. 



parameters 6RLM(t), describing the octupole surface vibrations, approximately 
satisfy an equation of motion of the damped-oscillator kind |24| : 

D 3 5R 3M (t) + l3 5R 3M (t) + C 3 SR 3M (t) = . (96) 

The coefficients D 3 , 73 and C3 are easily evaluated, showing that the oc- 
tupole surface oscillations described by Eq. (55)1 are overdamped. Another 
interesting result of this semiclassical analysis concerns the shape stability of 
heavy spherical nuclei against octupole-type deformation: within the present 
model the spherical shape is stable. 
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4 Conclusions 



The linearized Vlasov equation, that can be seen as a particular case of the 
Landau kinetic equation for the phase-space-density fluctuations, gives a good 
qualitative description of the low-energy isoscalar nuclear response of different 
multipolarities. This collisionless equation, which has been initially derived for 
other systems, can be applied also to nuclei because in nuclear matter the mean- 
free-path of nuclcons close to the Fermi surface is larger than the typical nuclear 
dimensions. This fact has two consequences: 

• the mean-field approximation is a reasonable one in the study of low- 
energy nuclear response, 

• finite size effects are important and should be taken into account. 

Hence, because of the interplay between nucleon mean-free-path and nuclear 
dimensions, the Vlasov equation can be used to study the nuclear response to 
a weak driving field of long wavelength. Clearly in different physical situations, 
like those realized in the collisions of heavy ions of intermediate or high energy, 
collisions become more important and should be taken into account. Thus 
in nuclei finite-size effects are more important at low energy, while collisions 
between nucleons become more and more important with increasing energy. 

As pointed out by Kirzhnitz and collaborators jxj, in finite systems, the 
boundary conditions satisfied by the fluctuations of the phase-space-density be- 
come essential. While in quantum mechanics these boundary conditions are au- 
thomatically enforced by the requirement that the wave function of bound states 
decreases exponentially outside the system, in the semiclassical kinetic-equation 
approach there is no similar requirement and the boundary conditions satisfied 
by the phase-space-density fluctuations must be imposed by using some reason- 
able criterion. Here we have studied the small fluctuations of the phase-space- 
density induced by applying a weak external driving force to spherical nuclei 
and have assumed a sharp-surface model for the spatial density and mean field 
in heavy nuclei. Then we have compared the experimental isoscalar strength 
functions with those calculated by imposing two different kinds of boundary 
conditions on the phase-space-density fluctuations: fixed- and moving-surface 
boundary conditions. On the whole we find a better qualitative agreement with 
experiment for the moving-surface response. 

As a final comment, we would like to add that the present semiclassical 
theory can be applied also to other systems of many fermions in which it is 
important to take into account the finite size. Atomic clusters and magnetically 
trapped droplets of fermions are interesting examples. Clearly, quantum cal- 
culations are more appropriate, however it might be of some interest to see to 
what an extent the quantum features stand on a classical "skeleton", and this 
can be appreciated most clearly within the present approach. 
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Appendix A: Moving-surface response function 

In this Appendix we give a few more details on the derivation of the analytical 
expression (|72(l for the multipole response function. 

The solution SfjJ^j of the linearized Vlasov equation with the boundary 
condition l|64l) can be written as [cf. Eq. (|29p] 



±*+"<r>»)[ f dr'B^ N (r')e^ N(r '^ + C± (e, A, w)] (97) 
J n 



with the functions B MN given by an equation similar to Q24JI 

B^ N (e,X,r,cj) = F'(e)[^- ± tN S ^]\PQL M (r) + SV£&(r,u)] (98) 
and the functions C± given by 

e 2i<t> N (R,uj) jj _ jj ^ 

(99) 

The mean- field fluctuation 5V™ M (r, u)) is a crucial quantity in our calcula- 
tions. Usually phenomenological models that describe the physical properties 
of the medium differ mainly in the assumptions made about this term. Our 
present approach is no exception to this general rule and Eq. (I72|) has been 
derived by assuming that 



5Vl n Mr,cj)=pK L r h 1l v L (cj), (100) 



with 



and 



ft£(w) = i / *//■/'" /*'' <} o i \i (ion 



de I d\ A 



u r (e,A,r) 



5/^+(e,A ) r ) w)+5/^ v (e,A,r,o;) . (102) 



The function (|101|l is not the response function l|72|l because, as discussed 
in |23| . in the moving-surface case the response function should include an ad- 
ditional term in order to take into account the shape changes, thus a more 
satisfactory definition of the multipole response function in the moving-surface 
case is (see also |3"2"] 1 

K l (lu) = i j drr 2 r L 8g LM (r,iu) , (103) 
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with 

Sg LM (r,uj) = Sg LM (r,oj) + g S(r - R)8R LM {w) , (104) 

giving 

il L {u) = nl(u) + ^R l+2 q q SR L m(lo) . (105) 

The response function (|72|l corresponds to l|105|l . rather than to l|101|) . The 
equilibrium density go appearing in Eq. (|104|l is go — -^(pp/K) 3 . 

In order to obtain the explicit expression H72[) of the response function 1)105(1 , 
we need the explicit expressions of the function R}^(lo) and of the collective co- 
ordinates 5Rlm{<+>). For deriving these quantities, the moving-surface solution 
ij^TI) should be expressed in terms of R}[{lu) and of SRlm{^)- By replacing the 
mean-field fluctuation (|100|l into the quantities B^ N and C±, given by Eqs. 
1(98(1 and (|99|l . and by inserting the resulting expressions into Eq. (|97(l . the 
fluctuation 5f^ N can be written as 

8f*f N (e,\,r,u) = */°f±(e,A,r,w)[l + 

e ±i* N (r,w) 

+ f, ( £ ) ■ r^, r5 „ rnv r (e,\,R)w6R LM {u) , (106) 

with the zero-order solution Sf®^ given by Eq. (|45|l . Now, by inserting the so- 
lution I|1U6|) into Eqs. (|101fl and l|71|l . we obtain a system of algebraic equations 
for the functions lZ\{u)) and 8Rlm{^>), that can be written as 

7#M = fc£(o;)[l + KL ^H]-ifx LH + eoi? i+3 l^#^ (107) 



7? 



and 



«?lmM - ^-{/3i? x £M[l + kz^H] (108) 
+3f£(w)«i2 £ wH + l3K L g a R L+i nY(cu)Y 

with the functions K° L {u), X°M and Xi( w ) S iven b y E q s - ISJ, (EJ and 
respectively. Solving the system (|107l - 1108(1 gives the explicit expressions of 
Ti)[(uo) and 8Rlm{^), which read 

^(o;)=^H-- ^— (109) 

and 

i^H „ xlH + K L g R L+3 K L (iv) 

i? p [c L -XLH][i- K LnH] + KL [x° L M + eoi? i+3 ] 2 ' 1 ! 

The fixed-surface collective response function 7Zl(oj) appearing in Eq. (|109|) is 
given by Eq. 1(53(1 . Finally, by inserting Eqs. 1(109(1 and 1110() into the response 
function (|105J) and taking into account that goR 3 = -j^A , we find Eq. (|72|l . 
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Appendix B: Fourier coefficients 



In this Appendix we collect the expressions of the integrals l|77|l needed here. 
Since the spherical harmonics in (|78|l vanish unless N has the same parity al 
L, we only need the corresponding integrals. Normally we need the coefficients 
with k = L, however for compression modes we also need k = L + 2. The 
coefficients involved in the monopolc response are 

L = 0, N = 

Q { nO +2 \x) = ^ [ R dr r \ cos# nA Kr)] (111) 
J J ri v r (eF,A,r) 

= {-) n R 2 „ 2 , . for n^O 

S n N ( X ) 

2 

= R 2 (l--x 2 ) for ra = 0, 

while for the quadrupole response they are: 

L = 2, N = 
same as monopole , 

and 

L = 2. N = ±2 



^' = <-'" fl2 ST«( 1 + w S- < 112 » 



zN (x) ^ ' s nN (x) 
In the dipole case we need the coefficients with 

L = l, N = ±l 



Q { nk*) = (-r R ^M ( 113 ) 



1 



for translation modes and 



- (-)"« 3 i(i + - -At) (H4) 



^jvfc) ^ ' 3 ^' S nN {x) S 2 nN (x) 



for compression modes. 

Finally for the octupole response we need: 

L = 3. W = ± i 5± 3 



4 .rVl 3 



:ir 



(115) 
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For a given nucleus, the integrals Q n h could depend on two variables: the 
nucleon energy cf and its angular momentum A. For the square- well potential 
however, they display a scaling property and depend only on the variable x. 
Moreover their ^-dependence factorizes because Q^n oc R k . As a consequence 
the A-dependence factorizes also in the zero-order propagator l|75|) . that takes 
the form of an A-dependent factor times a universal propagator. 
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